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Abstract 

We present the first model independent determination of S- and P-wave helicity amplitudes An, n = 0, 1, 
A = S,L,U from CERN measurements of Tr~p 7r~7r+n on polarized target at small t and dipion 
masses 580 -1080 MeV. The purely analytical determination of the helicity amplitudes is made possible 
by our finding analytical solutions for relative phase = ^SaU) ~ *J's„(*) between 5- wave amplitudes 
of opposite transversities u = up and d = down for each set of the solutions of transversity amplitudes 
Au{i), Ad{j),i, j — 1,2. Of the six possible solutions for Uij only the solution with ujij — tt yields physical 
helicity amplitudes Ao{ij), Ai{ij),i, j = 1,2. Assigning p°(770) phase to the dominant P-wave helicity 
flip amplitude Li{ij) necessitates a phase of the S'-wave helicity flip amplitude Si{ij) that is near to the 
p'~'{770) phase. Both amplitudes show resonant structures around 980 MeV for all solution sets i,j — 1,2. 
These amplitudes are thus consistent with p°(770) — /o(980) mixing observed previously in the reduced 
transversity amplitudes. The relative phases ujij — tt satisfy certain self-consistency condition that must 
be satisfied in order for the four sets of solutions Au{i), Ad{j),i, j — 1,2 to be all physical amplitudes 
that can be identified with co-evolution amplitudes describing the interaction of the pion creation process 
with a quantum environment. This test on the phases LUij provides a new test of Kraus representation of 
the mixed final state density matrix in ir^p — s- 7r~7r+n and connects it to the experimentally measured 
amplitudes. This connection validates further the view of pion creation process n^p 7r~7r"*"n as an open 
quantum system interacting with a quantum environment. The probabilities pij determining the final 
mixed state p/ = ^PijPfiij) in terms of solution states Pf{ij) are experimentally measurable in measure- 
ments of recoil hyperon polarization in self-analyzing processes 7r~p — > tt~K~^A'^ and K~p — > tt^tt"*" 
on polarized targets. The probabilities pij provide information about the quantum state of the environment. 
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I. INTRODUCTION. 



One of the fundamental assumptions of Quantum Theory is the unitary evolution of quantum 
systems. The assumption means that the quantum system S is isolated and its future state can 
be calculated. Quantum Theory admits non-unitary evolution of open quantum systems that are 
not isolated but interact with a quantum environment E. In this case it is the system S plus the 
environment E that together undergo a unitary coevolution [l|, Q, 0] • 

Theory of elementary particles is based on the assumption that the particle interaction is an 
isolated event in the Universe subject to unitary evolution [E, 0]- In 1982 Hawking questioned 
the universal validity of the unitarity assumption in particle interactions 0, (sl . He suggested that 
particle interactions are processes that interact with an environment of quantum fluctuations of 
space-time metric which induce a non-unitary evolution of the particle interaction - at any energy. 
As the result of the non-unitary evolution pure initial states evolve into mixed final states. In 
contrast, a unitary evolution evolves pure initial states into pure final states. 



Spin physics [9] is ideally suited to test Hawking's ideas. In our previous work [lu] we 
investigated the unitarity assumption in the pion creati on p rocess tt~p — > 7r~iT~^n measured at 
CERN on polarized target at 17.2 GeV/c [ll|, [l3, [l3, Q Following the initial work by Lutz 



and Rybicki [1,1] we used a spin formalism to relate final state density matrix to transversity 
amplitudes with definite dipion spin, helicity and naturality. Imposing the purity of the final 
states for specific initial pure states as required by the unitarity leads to a set of constraints on the 
transversity amplitudes. The unitarity constraints are violated by measured amplitudes at large 



momentum transfers t obtained in a model independent amplitude analysis of CERN data 15 1. 
This result suggests that in Tr~p vr^vr^n pure states evolve into mixed states and that the 
pion creation process behaves as an open quantum system interacting with a quantum environment. 

There are other indications of a non-unitary evolution in ir^p vr^vr+n. At any dipion mass 
m, the S- and P-wave subsystem of reduced density matrix measured at polarized target is ana- 
lytically solvable in terms of reduced S-and P-wave transversity amplitudes with definite t-channel 
naturality [lH |. There are two physical solutions Au{i),i = 1,2 and Aii{j),j = 1,2 for transversity 
amplitudes with target nucleon transversity "up" and "down", respectively Associated with 
each set of the amplitudes is a final state density matrix pf{ij). Each Pf{ij) carries specific infor- 
mation about the pion creation process and selecting one of the four sets of solutions amounts to 
a loss of such information. We can retain all the information if we assume that the measured final 
state density matrix is a mixed state of the solutions Pf{ij) 

P/ =PiiP/(ll) +Pl2/3/(12) +P21P/(21) +P22P/(22) (1.1) 

where the probabilities YlPij ~ this case the pure initial states will evolve into mixed final 

states even when Pf{ij) are pure states. 

Observation of p^{770) — /o(980) mixing in ■K~p — > 7r~7r"'"n provides another evidence for 
non-unitary dynamics of this process. Previous amplitude analyses of CERN data at small and 



large t found a rho-like resonance cr(770) in the 5- wave amplitudes 13|, llSl . llTI . llSl . Il9l |. In a new 
high resolution amplitude analysis of CERN data at small t [3| we show that this resonance 
must be p^{770) and that data indicate ^"(770) — /o(980) mixing also in the P-waves. This 
interpretation explains why no rho-like resonance has been observed in S'-wave amplitudes in 
TT~p — > vr'^vr'^n where P- waves do not contribute In Ref. [3] we propose to account for 

the ^"(770) — /o(980) mixing by introducing a new CPT violating interaction of the produced 
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resonant qq modes with interacting degrees of freedom of the environment which are identified 
with quantum numbers labehng the solutions for the amphtudes. 

In Section II. we show that the requirement (1.1) that ah four sets of solutions for transversity 
amplitudes 74ti(i), ^^(j), i, j, = 1,2 are accepted physical solutions leads to a self-consistency 
condition among relative phases LOij = ^s^ij) — *^5u(*)- In Section III. we introduce Kraus 
representation of the mixed final state density matrix (1.1) arising from a unitary co-evolution 
of the pion creation process with a quantum environment. We define co-evolution amplitudes 
which are identified with the measured solutions for transversity amplitudes. The co-evolution 
amplitudes must satisfy the self-consistency condition so that its experimental test is a test of 
validity of Kraus representation in pion creation process 7r~p 'K~'K^n. 



It has been generally believed that the determination of the relative phase oj requires difficult 
measurements of recoil nucleon polarization which involve the required bilinear terms between the 
amplitudes with nucleon transversities "up" and "down". In this work we show that the phases 
LVij can be determined from measurements on transversely polarized targets using relationships 
between helicity and transversity amplitudes. In Section IV. we present the helicity amplitudes 
An{ij),n = 0, 1 and the expressions for their moduli and bilinear terms in terms of known bilinear 
terms of measured reduced transversity amplitudes and unknown cos a; and sin a;. In Section V. 
we impose a self-consistency requirement between the product of two moduli squared and the 
square of magnitude of the corresponding bilinear term. This trigonometric condition yields 6 
solutions for the pairs cosw and sinw for each solution set i,j = 1,2. 

The results are presented in Section VI. .Two solutions with cos uj ^ and sinu; 7^ are 
rejected because the corresponding solutions for helicity amplitudes exhibit unphysical chaotic 
behaviours and lack a clear resonant structure at ^^^(770) mass in P-wave amplitudes. Three 
solutions with cosw = 0, sinu; = ±1 and cos a; = +1, sinu; = are also rejected because the 
corresponding helicity amplitudes do not satisfy the requirement of pion exchange dominance 
of helicity flip amplitudes. The remaining solution with coso; = — 1, or w = vr, yields a unique 
solution for helicity amplitudes that satisfy the requirement of pion exchange dominance and 
shows resonant structures in the P-wave. Assigning ^^^(770) phase to the dominant P-wave 
helicity flip amplitude Li{ij) necessitates a phase of the S'-wave helicity flip amplitude Si{ij) 
that is near to the p^{770) phase. Both amplitudes show resonant structures around 980 MeV 
for all i,j = 1,2. These amplitudes are thus consistent with p°(770) - /o(980) mixing observed 
previously in the reduced transversity amplitudes {l^. 



The physical solutions for the relative phases tOij = vr satisfy a self-consistency condition that 
must be satisfled by co-evolution amplitudes. The transversity amplitudes Au{i), Afi{j),i, j = 1,2 
can thus be identified with co-evolution amplitudes connecting the Kraus representation to the 
experimentally measured amplitudes. This connection validates further the view of pion creation 
process ir'p TT~iT'^n as an open quantum system interacting with a quantum environment [id ]. 



In Section VII. we show explicitely that the unique solution for iVij = vr leads to a unique solution 
for the complete set of S- and P-wave density matrix elements defining each state Pf{ij) which 
results in a unique form of Kraus representation of the mixed final state (1.1). In Section VIII. we 
show how measurements of recoil hyperon polarization in K~p — > 7r~7r"'"A'' and 7T~p ir^K^A^ 
on polarized target can be used to determine the probabilities pij in the mixed final state (1.1). 
The probabilities pij provide information about the quantum state of the environment. The paper 
closes with a summary in Section IX.. 
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II. S- AND P-WAVE AMPLITUDE ANALYSIS OF vr p ^ tt 7r+n ON TRANSVERSELY 
POLARIZED TARGET AND SELF-CONSISTENCY OF THE CENTRAL HYPOTHESIS. 

Amplitude analysis is an integral final part of any measurement of 7r~p — > TT~TT'^n on 
transversely polarized target when recoil nucleon polarization is not observed. It is a model 
independent conversion of the measured reduced density matrix elements into moduli and phases 
of nucleon transversity amplitudes with definite dipion spin J, helicity A and t-channel naturality 
V 



16l |. For any dipion mass m the S- and P-wave subsystem of the reduced density mantrix is 



analytically solvable in terms of reduced transversity aniplitudes A and A with target nucleon 



transversity t = u (up) and t = d (down), respectively [Ig]. In our notation jl6(] A = S,L are 
unnatural exchange S- and P-wave amplitudes with helicity A = while A = U,N are unatural 
and natural exchange P-wave amplitudes corresponding to combinations of A = ±1 [lo| . 



In a previous work [16(] we have shown that amplitude analysis of the complete 
S- and P-wave subsystem of the reduced density matrix yields four sets of solutions 
A{i), N{ij), A{j), N{ij),i, j = 1,2 where A = S,L,U. This requires measurements with 
target polarization not only with transverse component Py but also with planar components Px 
and Pj. Such measurements are feasible with modern polarized targets [91. 

For any set of solutions i, j = 1,2 the reduced transversity amplitudes A, A are defined in terms 



of transversity amplitudes Au{i),A(i{j) as follows 16|] 



Sii) = \Suii)\, SU) = mj)\ (2.1) 
L{i) = \Lu{i)\expi{^LAi) - <l>s^{i)) , L{j) = |Ld(j)| exp i ($L,(j) - ^5d(i)) 
U{i) = |C/,(z)|expi($j;„(z) -<I.s„(i)). U{j) = |i7rf(j)|expi($t;,(j) -^5,(j)) 

N{ij) = \N^(i)\ expi ($ivji) - ^sM)) , = |iVd(j)| exp (i$jv,(j) - %„(i)) 

where is the phase of the amplitude At-. The reduced transversity amplitudes are related to 
transversity amplitudes by phase factors 

A„(i) = A(i)expi$5ji), Ad(j) =]4(j)expiu;ijexpi$5„(i) (2.2) 
for unnatural exchange amplitudes A = S,L,U and 

Nuii) = N{ij)expiuJijexpi<^sA^), N^ij) = N{ij) exp i<^>s^ii) (2.3) 
for natural exchange amplitude N. In (2.2) and (2.3) ^^^(i) is the arbitrary absolute phase and 

'^ij = ^^S.U) - ^sAi) (2.4) 

is the relative phase between 5-wave amplitudes of opposite transversity. The moduli 
— \^u{i)\ and |A^(ij)| = \Nii{j) depend only on the indices i and j, respectvely, 
and are determined by measurements on transversely polarized target. In contrast, the phases 
of N{ij) and N{ij) depend on both indices and are determined by measurements with planar 
components of target polarization. 



5 



Up to now it has been generally believed that experimental determination of the relative phases 
LVij requires difficult measurements of planar components of recoil nucleon pola rization which 



ucieon pola r iza 

involve the interferences between transversity amplitudes Au and 11 , E^. lid. EtI . Il8l |. Each 
set i,j = 1,2 of reduced transversity amplitudes would determine the corresponding phase coij. 
In Section IV. of this work we shall show that the phases ujij can be determined analytically in a 
process of converting the reduced S- and P-wave transversity amplitudes A and A, A = S, L,U 
into S- and P-wave helicity amplitudes An with definite t-channel naturality where n = 0, 1 is 
nucleon helicity flip. This process also involves the necessary interferences between transversity 
amplitudes Au and A^. 

Our central hypothesis (1.1) that all sets of transversity amplitudes Au{i), Aii{j),i, j = 1,2 are 
valid physical solutions imposes experimentally testable self-consistency constraint on the relative 
phases Uij. Taken separately, each set Au{i),Ad{j) has its own absolute phase $s'„(^) and Uij are 
in general not related. Taken together, all sets Au(i), Ad{j),i, j = 1,2 can have only one absolute 
phase, say ^^^(l). Then we can write for A = S, L,U 

Au{l) = exp(icD5jl))^(l) (2.5) 

Au{2) = exp(icD5jl))exp(iO^(_2) 

Adil) = exp(i$5jl))exp(iwii)]4(l) = exp(i$5jl)) exp(iO exp(iw2i)^(l) 

Ad{2) = exp(i$5jl))exp(iwi2)^(2) = exp(i$5„(l))exp(iOexp(iw22)^(2) 



where 



e = $5j2) - cDsJl) (2.6) 



Prom the definitions lou = $5^(1) — $5^(1) and UJ21 = <I>5^(1) — <I>5'„(2) it follows that uJn —0J21 = C 
while from the definitions u;i2 = ^"5^(2) - <I>5„(1) and UJ22 = ^sA'^) - ^s„(2) it follows that 
a;i2 — u:22 = ^- The phases cOij are not independent but satisfy a condition 

^ = wii - UJ21 = u)\2 - UJ22 (2.7) 

For each set of reduced transversity amplitudes A{i)^A{j), A = S, L,U the conversion to helicity 
amplitudes provides an independent determination of the corresponding phase cuij for the set. If 
all the transversity amplitudes Au{i), Ad{j),i, j = 1,2 form a self-cosistent set of valid physical 
solutions then the condition (2.7) must be satisfied. Any violation of the condition (2.7) means 
that the calculation of uju and UJ21 require phases $5„(1) and $5^(2) different from those in 
the calculation of uji2 and uj22- In such a case the four solutions for the transversity amplitudes 
Au{i),Ad{j) do not form a self-consistent set of valid physical solutions. The condition (2.7) is 
thus an essential test of the central hypothesis (1.1). 



III. KRAUS REPRESENTATION AND ITS TEST IN tt p ^ tt 7r+n. 

In order to introduce the central concept of co-evolution amplitudes we ffist briefly review the 
Kraus representation for for reduced density matrices of open quantum systems interacting with 
an environment. It is the co-evolution amplitudes which involve the interacting degrees of freedom 
of the environment and with which we shall identify the solutions for the transversity amplitudes. 

The co-evolution of an open quantum system S with a quantum environment E' is a unitary 
evolution [l|, [3] 

Pf{S, E) = Upi{S, E)U+ = Upi{S) ® p^{E)U+ (3.1) 
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The initial state of the environment is in general a mixed state 

Pi{E) = ^Pu'\et X eit\ (3.2) 

where |e£ > are quantum states of interacting degrees of the environment and '^pu = 1- The 

Hilbert space of the environment has a finite dimension. It is given by a condition d\m.H{E) < 
dim Hi(S) dim Hf{S) [2l]. Here H(E), Hi(S) and Hf{S) are Hilbert spaces of the environment E 
and of the initial and final states of the system S. After the interaction the system S is fully 
described by reduced density matrix Pf{S) given by Kraus representation 

pj{S) = TrE{pf{S,E)) = '^'^PmnSlmPi{S)S^^ (3.3) 

where the operators Sirn =< ^i\U\em > satisfy a completness relation ^ ^ S^^Si^ = I- Kraus 

representation assumes that the initial states Pi{S) and pi{E) are separable to ensure complete 
positivity of Pf{S) fJ]. 

In our next step we associate the two solutions for transversity amplitudes Au{i) and vlrf(j), 
i,j = 1,2 with two single qubit states \i > and \j >, respectively. Then the hypohesis (1.1) 
allows us to identify the four degrees of freedom of the environment \e£ > allowed by the condition 
dimH{E) < dimHi{S) dimHf{S) = {2sp + l)(2s„ + 1) = 4 with the four two-qubit states \ei >= 
\i > \j >. The co-evolution amplitudes are then defined by matrix elements 

C//'^(^m) =< JXv,Tn\ < ee\U\em > |Orp > (3.4) 

where J, A and r] are dipion spin, helicity and t-channel naturality, and t = Tp and are transversi- 
ties of target proton and recoil neutron, respectively. Since the transversity amplitudes can possess 
only one solution at a time, the co-evolution amplitudes U^^^{tm) must be diagonal for any dipion 
spin J and naturality rj 

U^^lilm) = ui'^^{U)6,m = Ui'^^iij, ij)Sij,i,f = Uil{i3)5ij,,r (3.5) 

where 

Uiliij) =< J^, rn\ < ij\U\ij > \0u > (3.6) 



UxHiiij) =< J>^V,rn\ < ij\U\ij > \0d > 
The requirement that the Kraus representation leaves invariant the spin formalism used in the 



data analysis necessitates 10(] that the co-evolution amplitudes C/_^^(ij) transform under P-parity 
as a two-body P-parity conserving process vr" +|? — > " J(7r~7r"'")" + n with parity P = (— I)'' for 
the dipion states " J(7r~7r+)" [10]. This means that there is no vector associated with the quantum 
states \i > \j > of the environment. In particular, there is no energy-momentum associated with 
these quantum states. The interacting hadrons conserve their energy-momentum and there is no 
exchange of energy-momentum with the environment, in agreement with the original proposal by 
Hawking for particle processes interacting with quantum fluctuations of the space-time metric [7|]. 
Instead, the interaction with the environment is a non-dissipative dephasing process. The co- 
evolution amplitudes Uji can then be written in a form 

Ufi = If, + i{27,f6\Pf - Pi)Tfi (3.7) 
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where Pi and Pf are total four-momenta of the initial and final hadron states and Tfi is the 
transition matrix for the process \7r~p > +\i > \j >^ |7r~7r"*'n > +\i > \j >. 

In a final step we identify the co-evolution transtion amplitudes with solutions for transversity 
amplitudes of any dipion spin 

Txliij) =< J^r,, Tn\ < ij\T\ij > \0u >= Ailii) (3.8) 

Tx^iiij) =< J^V,rn\ < ij\T\ij > \0d >= Aj^(i) 

where A'^^_^{i), A'^^^{j), i,j = 1,2 are the experimental solutions for transversity amplitudes. The 

set of co-evolution amplitudes T^^{ij) is a consitent set with a single absolute phase <I>q~^(11) 
that satisfies the consistency conditions (2.7). The consistent identification of the co-evolution 
amplitudes with the experimentally determined solutions for transversity amplitudes requires that 
the latter also satisfy the conditions (2.7) so that <I>q^^(11) = <I>5^(1). Since we have identified 
the interacting degrees of the environment with the quantum numbers labelling the solutions for 
transversity amplitudes, the test of condition (2.7) is de facto a test of Kraus representation of 
the final state Pf{S) in pion creation process n^p — > Tr^vr^n. 

In our previous work [3] we have introduced a quantum number g = ±1 to label the solutions 
for the amplitudes. Instead of using the solution qubits \i > and \j > to define the quantum 
states \i > \j > of the environment and the co-evolution amplitudes T^^{ij), we could have used 
the qubits |g„ > and \gd > to define equivalent states \gu > \gd > of the environment and the 
equivalent co-evolution amplitudes T^^^{gugd) = A'!^^{gT). These states more closely reflect the 
qubit nature of the interacting degrees of the environment and may possess a deeper physical 
meaning. However, in the context of the present work we prefer to use the labels \i > and \j > for 
the solutions. 



IV. HELICITY AMPLITUDES AND THEIR BILINEAR TERMS. 

Helicity amplitudes with definite f-channel naturality were defined in Ref. l3, 11 1 for any dipion 
spin J and helicity A and their relations to transversity amplitudes of definite f-channel naturality 
were given. The S- and P-wave helicity nonflip and flip amplitudes and Ai are related to 
transversity amplitudes A,-, r = n, d, A = S^L,U,N hy relations 

A„ = i-^(^. + (-l)M,) (4.1) 

where n = 0, 1. In terms of reduced transversity amplitudes we can write for the unnatural 
exchange amplitudes A = S, L,U 

An = + i-irAexp{iu))exp{z<^sJ (4.2) 

where u; = $5^ - . 

Measurements on transversely polarized target yield information on reduced density matrix 
elements and Py. Amplitude analysis of this data results in four sets of solutions for the moduli 
of reduced transversity amplitudes |A(i)|, |^(j)|, i, j = 1,2. Each set comes with a fourfold sign 



8 



ambiguity in the phases of the ampUtudes {A{i), A{j))f^ where ee = ++, H — , — h, are the signs 

of relative phases $ls'(0 = ^L„(«) - ^5„(0 and ^lsU) = ^LaU) " ^5^0')- These phases do not 



change sign as a function of dimeson mass m in any of the 8 reactions analysed in Ref . 16(] . The 



change of sign of <^L5(i) {^Ls{j)) which is the phase af the amplitude L{i) {L{j)) results in the 



change of sign of the phase of amplitude U{i) {U{j)) [16t], or complex conjugation of amplitudes 



A{i) {A{j)). We can thus write the four sets of phases for each z,j 

= (^(i),^(i)) (4.3) 
= {A{i),A{jy) 

{A{i),A{j)).^ = {A{i)\A{j)) = {{A{i),A{j))^^r 
= {A{i)\A{jr) = {{A{i),A{j))++r 

The sets {A{i),A{j)) |_ and {A{i),A[j)) are complex conjugate of sets (^(i), and 

{A{i), A{j))j^j^, respectively. Only the measurements with planar target polarization can select 
uniquely one of the solutions for the phases [3l • 

For any given set of signs ee of the phases we can now calculate the corresponding helicity 
amplitudes An{ij)ee,'i', j = 1,2. Using the expressions (2.5) and omitting the indices ee for brevity, 
the helicity amplitudes read 

(4.4) 



Anill) = 


^/2 


+ (- 


-l)M(l)exp(ic^n))exp(i$5„(l)) 




An{l2) = 


^/2 


(^(1) + (- 


-l)"A(2)exp(ic^i2))exp(i$5jl)) 




Ani21) = 


V2 


(^(2) + (- 


-1)"A(1) exp(iw2i)) exp(i^) exp(i$5„ 


(1)) 


An{22) = 


{-iT 


(^(2) + (- 


-1)"A(2) exp(iu;22)) exp(iO exp(i'l'5„ 


(1)) 



It is easy to verify that the effect of change of signs of phases is again a complex conjugation of 
helicity amplitudes associated with a change of sign. Specifically, for any solution set i, j we obtain 

^o(ij)-+ = +(^o(u)+-)*, ASi)— = +{Ao{^j)++r (4.5) 

Ai{ij)-+ = -(^i(ii)+-)*, Ai(ij)__ = -{A,{ij)++r 
The change of signs of 5- wave phases results in change signs of relative phases lo and ^ 



Since with the phase sets and — + we do not get any new solutions, just complex conjugates 

of ++ and H — solutions, numerical calculations were done only for the sets with phases ++ and 
+-. 

Next we look at bilinear terms A^A*^ = |A„|2, A = S,L,U and AnB*,AB = LS, US, UL. For 
the sake of brevity we shall omit in the following the indices ij and ++, H — . Using (4.2) or (4.4) 
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we obtain 

I = 1 (l + \A\^ + {-\r2XA cos(a;) + (-l)"2y^ sin(a;)) (4.7) 

= \{Ia + (-1)"2Xacos(u;) + (-l)'^2yA sin(u;)) 

where Xa = Re{AA*), Ya = Im{AA*) and Ia = + is partial wave intensity. Note that 
Ys = Im{SS ) = as both 5 and S are real. For the bilinear terms AnB* we obtain 

AnB;^ = ^(aB* + AB* + {AB* e-'"^ + AB*e+''^)^ (4.8) 

The real part reads 

Re{AnB*) = ^ (Re{AB*) + Re(AB*) (4.9) 

{{Re{AB*) + Re(AB*)) cosw + {Im{AB*) - Im(AB*)) sinw)) 
The imaginary part reads 

Im{Ar,B*) = ^ (^Im{AB*) + Im(AB*) (4.10) 

+{-iy'{{Im{AB*) + /m(:45*)) cos a; - {Re{AB*) - Re(AB*)) sinw)) 

It is apparent from (4.2) and (4.4) that the knowledge of u allows to determine the helicity 
amplitudes up to an absolute phase for eaeh solution set i,j and phase set ++ and H — . As we 
show in the next Section, the phase uj can be determined analytically from the consistency condition 

\An\^\Bn\^ = {Re{AB*)f + {Im{AB*)f (4.11) 

where the terms are given by (4. 7), (4. 9) and (4.10). 



V. ANALYTICAL SOLUTIONS OF THE RELATIVE PHASE uj. 

In order to make use of the consistency condition (4.11) to determine oj , we first rewrite (4.7), 
(4.9) and (4.10) in a simplified form. To this end we define 

X{AB) = Re{AB*) + Re(AB*) (5.1) 

Y{AB) = Im{AB*) + Im(AB*) 
X{AB)+ = Re{AB*) + Re(AB*) 
Y{AB)+ = Im{AB*) + Im(AB*) 

X{AB)^ = ReiAB*) ~ Re(AB*) 
Y{AB)_ = Im{AB*) - Im(AB*) 



Then 



Re{AnB*) = ^(^X{AB) + {-l)''{X{AB)+cosLo + Y{AB)_smcj)^ (5.2) 
Im{A^B*) = ^(y (AB) + (-1)'' {Y{AB)+ cos CO -X{AB)_ sin Lv)^ 
Re{AnAl) = \Anf = ^(X{AA) + {-1^ {X{A A) + cos u + Y{A A) ^ sin u)^ 
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Note that X{AA) = Ia, X{AA)+ = 2Xa and Y{AA)^ = 2Ya- We can write 

1 
2 



ReiA^B:) = l-[XiAB) + {-irXGiAB)] (5.3) 



where 

Next we require that 
The l.h.s. of (5.5) reads 



Im{AnB*) = ^(Y{AB) + {-irYG{AB)) 
\Anf = ^[x{AA) + {-irXG{AA)') 
\Bn\^ = i(x(i?5) + (-l)"XG(i?5)) 

XG(AB) = X{AB)+cosuj + Y{AB)^smuj (5.4) 
YG{AB) = Y{AB)+cosu- X{AB)^ sin uj 

\Anf\Bnf = {Re{AB*)f + {Im{AB*)f (5.5) 



X{AA)X{BB) + XG{AA)XG{BB) (5.6) 

+{-l)'^X{AA)XG{A'A)X + {-l)'^X{BB)XG{A'A) 
The r.h.s. of (5.5) reads 

X{ABf + XG{ABf + Y{ABf + YG{ABf (5.7) 

+{-\Y2X{AB)XG{AB) + {-IY2Y{AB)YG{AB) 

Subtracting (5.5) with n = 1 from (5.5) with n = and using (5.4) we obtain equation hnear in 
cos a; and sin a; 

smu(x{AA)Y{BB)^ + X{BB)Y{AA)^ - 2X{AB)Y{AB)^ + 2Y{AB)X{AB)^^ = (5.8) 

- cos ui(^X{AA)X{BB)+ + X{BB)X{A'A)+ - 2X{AB)X{AB)+ - 2Y{AB)Y{AB)+^ 
which can be cast in the form 

sin a;W2 = — cos a; W"i (5.9) 

Using sin^ u + cos^ a; = 1 we find 

coso; = —r— (5.10) 
W 

where W = ^/wfTW^. 

Using (4.7), (4.9) and (4.10) it is straightforward to verify that the consistency condition (5.5) 
reduces to identity in the following two cases 

cosuj = 0, sinoj = ±1 (5-11) 
sinoj = 0, cosoj = ±1 (5-12) 

The first case leads to two solutions to = +7r/2 and lo = — 7r/2. The second case leads to another 
two solutions a; = and to = w. 
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VI. SOLUTIONS FOR S- AND P-WAVE HELICITY AMPLITUDES. 



A. Numerical calculations and their checks. 



Monte Carlo amplitude analysis of CERN data on Tr~p — > 7r~TT~^n measured on transversely 
polarized target at 17.2 GeV/c at small momentum transfers 0.005 < —t < 0.2 (GeV/c)^ was 
performed using a computer code AACERNM for dipion masses 580 -1080 MeV jl6l |. The mean 
values of moduli and phases of reduced transversity amplitudes were used as an input in a 
computer code HACERNMa to calculate uj and helicity amplitudes An,n = 0,1,^4 = S,L,U. 
The mean values of moduli and phases of the reduced transversity amplitudes satisfy strict 
normalization and phase conditions, res pec tively, and thus represent the true measured amplitudes 



at the i-bin average of 0.067 (GeV/c)^ [l6|,ll9(]. 



Equations (5.10) were used with amplitudes A = L, B = S to calculate uj{ij) for each set ++ 
and H — of signs of phases. Moduli lyln^J-B^p and interference terms Re{AnB*), Im{AnB*) for 
nonflip n = and flip n = 1 pairs of amplitudes AB = LS, US,UL were calculated to determine 
cosines and sines of the relative phases 

^Ar,B*) = <^>A„-'^>B„ (6.1) 

Since the calculations of uj, moduli and the interference terms are all entirely independent, the 
selfconsistency of the data and calculations was checked using the following three tests on the 
relative phases. The first test are trigonometric identities for n = 0, 1 for each pair of amplitudes 
AB = LS, US, UL 

cos^ ^{AnB*) + sin^ HAnB*) = 1 (6.2) 
The second test are phase conditions for n = 0, 1 

(^L„ - ^sj - (^a. - ^sj + (^{/„ - ^lJ = (6.3) 

The third test are cosine conditions for n = 0, 1 equivalent to phase conditions 

cos2^>(L„5;;)+cos2^>(;7„S*) + cos2$(C/„L;)-2cos$(L„5;)cos$(?7„5;)cos$(f/„L;) = 1 (6.4) 

All conditions are satisfied identically within the single precision calculation used by HACERNMa 
for all combinations of solutions i,j = 1,2 for both sets ++ and H — of signs of phases. For 
instance, typical deviation from 1 of the trigonometric identities (6.2) are of order 10~^ — 10~^. 



Errors on uj and helicity amplitudes were not calculated. This could be accomplished by combin- 
ing the core of AACERNM with the code HACERNMa and calculating the helicity amplitudes for 
each Monte Carlo sampling of the data that yields physical solutions for the reduced transversity 
amplitudes. The resulting distributions of values of helicity amplitudes would define their mean 
values, and the range of distributions would estimate their errors. These calculations are however 
beyond the scope and purpose of this paper. 



B. Unphysical solutions with cosoj ^ and sinuj ^ 0. 



The calculated values of uj{ij)-^-^ and uj{ij)^ show random variations as a function of dipion 

mass for all solution sets i,j = 1, 2 and signs of phases. To test the consistency condition (2.7) 

^ = UJu — ^21 = <jJi2 — ^22 (6.5) 
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we calculated the quantity 

A = (wii - UJ21) - (wi2 - W22) (6.6) 

The results are shown in Fig. 1. The large random variations of A indicate a strong violation of 
the consistency condition which requires A = 0. These results are thus inconsistent with Kraus 
representation and the assumption that all four sets of reduced transversity amplitudes are physical. 

One could hope that one of the solutions could be selected as a valid physical solution. How- 
ever, all amplitudes exhibit the same random behaviour and do not show the required resonant 
Breit-Wigner behaviour at p°(770) resonance. In Fig. 2 and Fig. 3 we present the results for the 
single flip amplitude Li for the signs of phases +- 1- and H — , respectively, with coscj = —W2/W 
to illustrate this behaviour. Moreover, the non-flip amplitude Lq has magnitude comparable 
to or larger than the flip amplitude Li, in contradiction with the physical requirement that Li 
dominates on account of pion exchange dominance at small t. The results with cos a; = +W2/W 
are similar. To ensure the positivity of ImLi the absolute phases $5„(^) were set equal to — tt. 

On the basis of the unphysical behaviour of the helicity amplitudes, all solutions with cos to = 
^W2/W are considered unphysical. They are excluded not only as co-evolution amplitudes but 
also as possible standard amplitudes. 

C. Unphysical solutions with cos a; = and sin a; = ±1. 

The two values for sinw = ±1 correspond to = ±7r/2 and Wjj^-i = ±7r/2 for all solution 

sets i,j = 1, 2 and signs of phases. Both solutions satisly the consistency condition (2.7). 

At small t the pion creation process 7r~p 7r~7r+n is dominated by pion exchange in the t 
channel. The pion exchange amplidues are the unnatural exchange helicity flip amplitudes »Si,-Li 
and Ui which must dominate the helicity non-flip amplitudes So,Lo and Uq, respectively. The 
CERN data correspond to small t bin average -t=0.067 (GeV/c)2. As the result, pion exchange 
dominance of the flip helicity amplitudes must be observed by the solutions. 

The solutions with sin a; = are excluded because they require that the non-flip amplitude 
Lq is larger than the flip amplitude Li. Because cosw = and the amplitudes Sq and 5*1 do not 
depend on sin a;, both solutions sinw = ±1 give |S'oP = l-S'ip = ^5/2 where Is is S'-wave intensity. 
The equal magnitudes of the 5-wave helicity flip and non-flip amplitudes contradict the pion 
exchange dominance for both solutions for u. As a result, both solutions are unphysical for any 
solution set i,j = 1,2 and signs of phases -|— |- and H — and are excluded from consideration as 
co-evolution or standard amplitudes. 

For the sake of comparison we present some results for sin a; = —1. Fig. 4 shows the amplitudes 
Si which are the same for both signs of phases -|— |- and H — . In Fig. 5 and Fig. 6 we show the 
amplitudes Li for the signs of phases ++ and H — , respectively. To ensure the positivity of ImLi 
the absolute phases $5u(i) were set equal to +7r/2 (— 7r/2 for sino; = +1). The imaginary parts of 
Li exhibit a clear /9°(770) resonant structure. There is no evidence of a chaotic behaviour in any 
of the amplitudes. 
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D. Physical solutions with coscj = 1 and sincj = 0. 



The two values for cosu; = ±1 correspond to <j-'jj^++ = 0,7r and u!ij^-\ = 0, vr for all solution 

sets i,j = 1,2 and signs of phases. Both solutions satisly the consistency condition (2.7) with 
^ = 0. From the definition (2.6) of ^ = <I>5^(2) — $5^(1) it follows that the amplitudes S„(l) and 
S'„(2) must have the same phase. 

The solutions with cosu; = +1 {lu = 0) require that < |5o|^ and < \Lq\'^ for 

any solution set i,j = 1,2 and signs of phases. The magnitudes of Si are small compared to 
magnitudes of So, in contradiction with the pion exchange dominance. These solutions are thus 
excluded as unphysical solutions. 

In the solutions with cos iv = —1 the non-flip and flip amplitudes are interchenged and the pion 
exchange dominance is observed. To ensure the positivity of ImLi the absolute phases ^^^(i) were 
set equal to — vr. The resulting solutions represent a unique solution for both helicity and transver- 
sity amplitudes (up to signs of phases to be resolved by measurements with planar components 
of target polarization). The amplitudes are consistent with Kraus representation and with the 
the central hypothesis (1.1) and can thus be identified with a unique set of co-evolution amplitudes. 

The results for helicity flip amplitudes are presented in Figures 7-11. The Fig. 7 shows 
the amplitude Si for both sets of signs -|— |- and H — of phases. While there is an indication 
for a p^{770) structure in solutions Si(ll) and Si(21), the solutions Si(12) and Si(22) show 
a broad structure indicating apparent suppression of p^{770). All solutions for Si show ex- 
pected resonant structures corresponding to /o(980) resonance. The Figures 8 and 9 show 
the amplitude Li for signs of phases -|--|- and H — , respectively, while the Figures 10 and 
11 show the same for amplitudes Ui. Both ImL\ and IrnUi show clear resonant peaks at 
p^{770) mass but have opposite signs. Their real parts are relatively structureless and small. 
Important structures appear in ImLi at /o(980) mass consistent with (770) — /o (980) mixing [ld |. 




The results for helicity non-flip amplitudes are presented in Figures 12-16. The Fig. 12 shows 
the amplitude So for both sets of signs -|--|- and H — of phases. The non-flip amplitude So is 
structureless and very small compared to flip amplitude Si. The Figures 13 and 14 show the 
amplitude Lq for signs of phases -|--|- and +— , respectively, while the Figures 15 and 16 show 
the same for amplitudes Uq. Solutions for both amplitudes appear structureless at p^{775) mass 
range with IitlLq changing sign around /o(980) mass. The non-flip amplitudes Uq are very small 
compared to flip amplitudes Ui. 

The solutions with cosw = — 1 represent the first model independent determination of all un- 
natural exchange helicity amplitudes in ■K~p — > 7r~7r"'"n. Assuming the resolution of sign ambiguity 
of phases in measurements with planar components of target polarization, the solution is unique 
without ambiguity. Moreover, assuming continuity of w and its constant value of vr at larger mo- 
mentum transfers t or dipion masses m, this solution can be selected in all amplitude analyses to 
provide a unique solution of co-evolution amplitudes at any t and m. 

E. The phases of amplitudes Li and 5*1. 

In the physical solutions with cosw = —1 the phase of amplitude Si is ^sA^j) = 90°. The 
phase of amplitude Li is shown in Fig. 17. It is nearly constant at ^ 90° below 950 MeV, which 
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appears suprising as one would expect its phase near the phase of p°(770) resonamce. Its sudden 
change of sign above 950 MeV is due to sudden change of sign in density matrix elements {p^)qs, 
{Pu)is, {Py)os, and {pl)is above this mass. 

We can assign to the amplitude Li the Breit-Wigner phase ^{p^) of p^{770) resonance by 
modifying the absolute phase — vr — vr + ^{p^) — ^Li{ij)- This uniformly modifies the phases of 
all amplitudes. In particular, the amplitudes Si acquire a modified phase 

(ij)' = (ij) + - (u) (6.7) 

Since ^Sii^j) — ^Li(u) ~ 0; the modified phase of Si becomes a somewhat modified phase of 
The modified phase of Si is shown in Fig. 18. The modification depends on the sign of the 
phases ++ and H — . 



For the sake of comparison we show in Fig. 19 the phase of Li for the solutions sinw = —1. The 
phases are again approximately constant below 950 MeV and can be assigned the Breit-Wigner 
phase of p'^(770). The modified phases of amplitudes Si are again near the phase ^{p^), as shown 
in Fig. 20. 



In 1997 Kaminski, Lesniak and Rybicki published a model dependent determination of helicity 
amplitudes based on CERN-Munich-Cracow analysis of the CERN data which did not determine 
the phases of the reduced transversity amplitudes but only their cosines [2l[ . They assumed Breit- 
Wigner phase for the amplitudes L and L and that the phases ^ls and ^ls change sign near 
p'^(770) mass. For the phases of amplitudes Si they obtained two "steep" solution and two "flat" 
solutions. Their "steep" solutions are near ^{p^) phase and are similar to our results in Fig. 18. 
Their "flat" solutions have much smaller rise and correspond closely to the CERN-Munich tttt 
phases shifts 5q. The authors emphasize in their paper that they can get the "flat" solutions only 
provided that the pha ses $£,5 and ^ls change sign near p^{770) mass. In our high definition 
amplitude analysis [l6|] we have determined the phases ^ls and and find no evidence for such 
a change of sign. Moreover, we show that all eigenvalues of the reduced density matrix are non-zero 
except at the point where the phase ^ls or change sign and where some of the eigenvalues 
vanish. Furthermore, the change of sign of the phases ^ls and causes abrupt and discontinous 
change of sign in all other phases [lo]. It is the unchangig sign of the phases ^ls and ^ls which 
prevents such anomalous behaviours of the eigenvalues and other phases which are not considered 
in Ref . 2l|] . We suggest that their analysis without the assumption of the change of sign of these 
phases would be similar to our model independent results. 



VII. UNIQUENESS OF KRAUS REPRESENTATION OF THE MIXED FINAL STATE 
DENSITY MATRIX. 

The unique solution for the phase iv = n implies a unique solution for the S- and P-wave 
transversity and helicity amplitudes from measurements on polarized target. In this Section we 
show that the four solutions for the S- and P-wave final state density matrices Pf{ij),i,j = 1,2 
are uniquely determined and, as a result, so is the Kraus representation of the mixed final state 
density matrix 

Pf =PllPf (11) +Pl2Pfil2) +p21/0/(21) +P22/0/(22) (7.1) 

Extension to higher dimeson masses where D-waves contribute will be given elsewhere. 
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For the specific case of vr p ^ tt TT~^n and similar processes with target polarization P = 
{Px,Py,Pz) the Kraus representation (3.3) takes the more explicit form 

Pf{e<p, p) = PUP fill, 6<p, P) + pi2P/(i2, e(t>, p) + P2iPf{2i, eel), P) + P22PA22, e<j), P) (7.2) 

where 

Pf{ij, e(p, P) = T{ij, e(p)pp{P)T{ij.e^)+ (7.3) 

In (7.3) T(ij,6(f)) is the matrix of co-evolution transition amplitudes defined by (3.7) 

<e(P,x\<ij\T\ij >\0,u> (7.4) 

where Vt = 9(j) describes the direction of vr^ in the dipion center-of-mass system, x is the recoil 
nucleon helicity, v is the target nuleon helicity and stands for incident pion helicity. The momenta 
of the particles are suppressed. In (7.3) Pp{P) is the density matrix of target nucleons 

Pp{P) = l{au + Pa) (7.5) 
where a are Pauli matrices and 0"^ = 1 is a unit matrix. 

Each final state density matrix pf(ij,6(j), P) is a single qubit density matrix corresponding to 



spin i of the recoil nucleon. It can be written in the form [2, [lO] 



Pfiij, 9^, P) = i(lO(ij, Oc^, P)a^ + e^, P)a) (7.6) 

where the traces I^{ij,d(f>, P) = Tr{a^ pf{ij,9(p,P)),t = 0,1,2,3 represent measurable intensities 
of angular distributions. Using pp = |(1 + Pa) we can write matrix elements of pf{ij,94>,P) in 
terms of components of target polarization 

{pf{ij,9cP,P))ll, = (p„(ii,e<A))H +P,.(p,(ii,^0))Jj +Py{py{zj,9cpm +PMij,OcP))l\ (7.7) 

In (7.7) the subscript u stands for unpolarized target P = 0. Using the decomposition (7.7) of 
Pf(ij,9(f),P) we find a decomposition for the intensities 

P{ij, 9(p, P) = Tr{a' Pfiij, 9^, P)) = 94>) + Pji{ij, 94>) + PyPy{ij, 94>) + Pji{ij, 9^) (7.8) 

where the components ll{ij,9(l)), £ = 0,1,2,3 and k = u,x,y,z of the intensities I^{9(t),P) are 
given by traces 

Ii{ij,9ct>) = Tr^,x"((<^')xx"(Pfe(^i,^</'))J4) (7.9) 

In general, a plane wave helicity state of two particles with helicities pi,p2 can be expanded in 
terms of angular helicity states 



\p9ct);prp2 >=Y,J^^^Dijcl),9,-cP)\pJX;pip2 > (7.10) 



where p is the momentum in center-of-mass system and J and A are the two-particle spin and 
helicity, and p = pi — p2- For two pions pi = P2 = and D'Iq{(I), 9, — </>) = Y/47r/(2J -|- l)Y^*{9, (p). 
The final state can be expanded in spherical harmonics 

\9(t>,x>=^Y^J*{9,(t>)\JX,x> (7.11) 

JA 
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TABLE L Density matrix elements Rep'^dPa / dtdm, Repyd^a/dtdm and Rep^d^a/dtdm, Repyd^a/dtdm in 
terms of reduced transversity amplitudes. Here aA — + and ta = \A\'^ — \A\'^ for A — U,N. The 
expressions are valid for any lu. 



pH' 


Repl 


pTs 
Pll 


\s\' + \s\^ 

\L\' + \L\' 


Pll 




Pi-i 


— ^{'^u — Cat) 




Re{LS* +LS* 



Repy RePu Rep. 



2 



\s? - \s\^ M)Ts -ip'u)Ts 

\L\'-\L\' -{P'y)hh -ipDhh 

U^U+Tn) -ip°y)ll+TN -{plVA+aN 

-^n/ - rjv) -ip°y)\li + TN -{pDHi + (TN 

Re{LS*-LS*) -ReipDl'i -Re{pl)ll 

V2Repl° Re{US* + US*) Re{US* - US*) -V2Re{pl)Yi -\/2i?e(p°)}° 

V2Repl\ RejLU* + LU*) RejLU* - LU*) -V2Re{p°y)ll -V2Re{pl 



TABLE IL Density matrix elements Impld^a/dtdm, Imp'^d^a/dtdm and Impl.d'^a/dtdm, Impj,d^ a / dtdm 
in terms of reduced transversity amplitudes. The expressions are valid for any uj. 



pii' 


Impl 


Imp^^ 




Impl 


Impl 




V2Impfi 
V2Impl\ 


Re{-SN* + NS*) Im{^ 
Re{-LN* + NL*) Im{~\ 
Re{+UN* - NU*) Im{- 


'SN* - 
'LN* ~ 
'UN* + 


NS*) 
NL*) 
-NS*) 


Re{SN* + NS*) Im{- 
Re{LN* + NL*) Im{~ 
Re{-UN* - NU*) /m(H 


-SN* - 
~LN* ~ 
-UN* + 


NS*) 
NL*) 
NS*) 



where J and A are dipion spin and helicity, respectively. Using (7.11) we obtain the angular 
expansion of transition amplitudes (7.4). The final expressions for if-iii, 9(1)) with independent 
angular density matrix elements assume P parity conservation and read [lOt] 

^^^'^■'^"^^ = ^ ^ 5;5;6J'CA(i^epi(^J))^^'i^e(y/(^0)yx^*(0,^)) (7.12) 

J< J' A>0 A' 

for {k, £) = {u, 0), {y, 0), (n, 2), {y, 2), (x, 1), {z, 1), {x, 3), {z, 3) and 

= ^ E E E CjMimpimii'iHyx'mY^:*m) (7.13) 

J< J' A>0 A' 

for {k,i) = (x,0),(z,0),(x,2),(z,2),(^z,l),(y,l),(tx,3),(y,3). In (7.12) and (7.13) Co = 1 and 
= 2 for A > 0. The factor Oj' = 1 for J = J' and CjJ' = 2 for J < J'. 



TABLE III: Density matrix elements Rep^d^a /dtdm, ReplcPa/dtdm and Rep'^d^a/dtdm, Rep^d^a /dtdm 
in terms of reduced transversity amplitudes. The expressions are valid for uj — tt. 



pjx ^^£Pi 

2Re{SS^) -2ImiSS_*) +(pi)™ -(pi)™ 
pll 2Re{LL*)_ -2Im{LL*)_ +{pl)}>h _ -(.P'^Voh 
pll Re{+UU* - NN*) Im{-UU* + NN*) +{pl)ll - 2Im{NN*) -{pDH - 2Re{NN*) 
pV-i Re{-UU* - NN*) Im{+UU* + NN*) +{pl)i-i ~ 2Im{NN*) -(pi)}Li - 2Re{NN*) 
Repll Re{L-S_*+Sr) Iwi-LT + ST) +Re{pl)ll -Re{pl)l° 
^Repll Re{US*+SU*) Im{-US* + SU*) +^/2Re{pl)ll -^/2Re{pl)Yl 
V^Repll Re{UL*+LU*) Im{-UL* + LU*) +V2Re{pl)ll -^/2Re{pl)ll 
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TABLE IV; Density matrix elements Imp^dPa / dtdm, Impyd^a / dtdm and Imp'^d^a/dtdm, Impyd^ a / dtdm 
in terms of reduced transversity amplitudes. The expression are valid for lu = n. 



pii' 


Impl 


Impl 




Impl 




Impy 




V2Imp°,l 


Re{SN* - 


'SN*) Re{SN* + 


'SN*) Im{^ 


-SN* ~ 


JN*) Im{- 


-SN* + 


SN*) 


V2Impl\ 


Re{LN* - 


LN*) Re{LN* + 


LN*) Im{- 


-LN* - 


LN*) Im{- 


-LN* + 


LN*) 


Imp]}ii 


Re{-UN* ^ 


-^UN*) Re{-UN^ - 


-UN*) Im{^ 


-UN* + 


UN*) Im{^ 


-UN^ - 


UN*) 



Lutz and Rybicki tabulated expressions for density matrix elements {Rl.)^^ = jtdm iPk)x\' 
terms of helicity amplitudes, helicity amplitudes with definite i-channel naturality, transversity 
amplitudes and transversity amplitudes with definite t-channel naturality [lH. Their tables are 
reproduced in the Appendix of Ref. ■ 

In Tables I.-IV. we present expressions for the complete set of measurable S- and P-wave density 
matrix elements pl.{ij) in terms of reduced transversity amplitudes A and A, A = S,L, U,N. The 
solution indices are suppressed for the sake of brevity. The expressions in the Tables I. and II. 
are valid for any phase uj and are not specific to the unique physical solution with uj = tt. The 
expressions in the Tables III. and IV. are valid only for the physical solution with to = tt. The 
expressions yield a unique solution for final state density matrices pf{ij,6(f), P). The mixed final 
state density matrix pf{6(j),P) given by the Kraus representation (7.2) is thus unique as well. We 
conclude that Kraus representation of the mixed final state is uniquely determined by measurements 
on polarized target alone. 



VIII. EXPERIMENTAL DETERMINATION OF PROBABILITIES 



A complete determination of the mixed final state pf{6(j),P) requires knowledge of the proba- 
bilities pij which are the diagonal elements of the initial state density matrix of the environment 

Pi{E) = J2 ^ Pij,i'j'\ij >< i'f\ (8-1) 

ij=l,2i'j'=l,2 

In this Section we indicate how the probabilities pij = Pij.ij can be determined in measurements 
of recoil hyperon polarization in processes such as K~p 7r~7r"^A^ and 7r~p it — K^hP using 
their self- analyzing weak decays. 

Recoil nucleon (hyperon) polarization vector Q{ij, 94>, P) is defined using a relation [lo| 

Qiij, 9^, P)l\ij, Q4>, P) = liij, 9^, P) (8.2) 

We can write 

Pf{ij, 9(P, P) = i(l + Q{ij, 6(1), P)a)l\ij, 9ct), P) = Pn{Q(i3))l\i3, ^'A, P) (8.3) 

The normalized final state density matrix p'j{ij,9(j), P) = pf{ij,9(j),P)/I^{ij,9(f),P) is simply the 
spin density matrix of the recoil nucleon PniQiij))- It will represent a pure final state if and only 
if the recoil nucleon polarization vector Q = {Q^,Q^,Q^) satisfies the condition IQp = 1 for all 
solid angles = (^,0) at any given dipion mass m and momentum transfer t 0, [s]. 
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Because of the linearity of the mixed final state density matrix (7.2), it has a form similar to 
(7.6) 

p;(00, P) = ^{l\9cP, P)a'' + m, P)a) (8.4) 

where 

l\e4>,P)= m3l'iij,0<P,P) (8.5) 

for ^ = 0, 1, 2, 3. The measured recoil nucleon polarization is similarly defined 

Q(00, P)/°(e0, P) = T{e<t), P) (8.6) 
so that we can write the mixed state in the form 

Pf{d<p, P) = i(l + Q{94>, P)a)l\9<P, P) = Pn{Q)l\e<p, P) (8.7) 
Inserting (8.3) into (7.2) and comparing with (8.7) we obtain 

Q(0<^, P) = V ( PmI^M:!^) Q^ij^ ecp, P) (8.8) 

In the dipion mass region where S- and P waves dominate the intensities Ocp, P) = l^{0(j), P) 

are the same for all i,j = 1, 2 [16] so that the measured recoil polarization is the average 

Q(^0,P)= ^ K,g(ij,0(/.,P) (8.9) 

Measurements on polarized targets determine uniquely the four sets of S- and P-wave density 
matrix elements and thus the four solutions of recoil polarizations Q{ij, 9(p, P) for any bin 
{9(1), m, t). Measurements of the recoil polarization Q{9(j), P) in four distinct bins {9(p, m, t) provide 
four independent linear equations for probabilities pij for each component of the polarization 
vector. The solutions for probabilities from all components must be consistent. The measurement 
of the probabilities pij provides information about the diagonal elements of the density matrix of 
the environment Pi{E) given by (8.1). 

Recoil nucleon polarization is not experimentally accessible in 7r~p — > 7r~7r"'"n. However, recoil 
hyperon polarization is accessible in measurements of K~p — > 7r~7r"'"A'' or 7r~p — > A9 through 

weak decays A'' — > p-K~ in measurements on unpolarized or polarized targets. Using the spin 
formalism for two-particle decays of a particle with spin ^ developed in Sections 8.2.1.(i) and 
8.2.1.(v) of the monograph on spin physics by E. Leader it is straightforward to derive a 
general form of normalized angular distribution W{9p, (f)p) of the protons in decays. It reads 

W{9p,(f)p) = — (\ + aQ^s\n9pCOS(pp + aQ'^si\\9pSm.(j)p + aQ^cos9r^ (8.10) 

where 9p, (j)p is the direction of the decay proton in center-of-mass system, a is a weak decay 
parameter measuring the real part of the interference between parity conserving and parity 
violating components of the decay amplitude and it is well known experimentally 

With sufficiently high statistics in {9(f), m, t) bins there will be sufficient statistics of A*^ decays 
to determine the recoil polarization vector Q{9(f>, P) with high precission. Simultaneous amplitude 
analysis of K~p — > 7r~7r^ AO or TT p ^ TT K~^A^ on polarized target will determine the four 
solutions for the recoil polarizations Q{ij,9(p, P), presuming the phase a; = vr is still the physical 
solution. In principle, these self-analyzing processes can thus determine the probabilities pij. 
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IX. CONCLUSIONS. 



We have presented the first model independent determination of unnatural exchange S- and 
P-wave helicity amplitudes in Tr~p TT~iT~^n from CERN measurements on transversely polarized 
target at 17.2 GeV/c for dipion masses in the range 580-1080 MeV and momentum transfers 
0.005 < —t < 0.20 (GeV/c)^. The analytical determination of helicity amplitudes is made possible 
by our finding of analytical solutions for the relative phase coij = ^SaU) ~ from the four 

sets of reduced transversity amplitudes A{i), A{j),i, j = 1,2, A = S, L,U found also analytically 



from the measured data [161 ]. 



The solutions for helicity amplitudes corresponding to two solutions for to with cosw ^ and 
smu> 7^ show chaotic behaviour and lack a clear resonant structure at p^{770) mass, and are 
rejected. Three solutions for helicity amplitudes corresponding to solution for u with sin a; = ±1 
and cos a; = +1 are also rejected since they do not satisfy the requirement of pion exchange 
dominance of helicity fiip amplitudes. This leaves a unique physical solution for helicity amplitudes 
corresponding to cos a; = —1, or tOij = vr. The solution is unique up to the sign ambiguity of 
phases of reduced transversity amplitudes to be resolved unambigously by measurements with 
planar target polarization [l^ . 



Assigning p'^(770) phase to the dominant P-wave helicity flip amplitude Li{ij) necessitates 
a phase of the S'-wave helicity flip amplitude Si{ij) that is near to the p°(770) phase. Both 
amplitudes show resonant structures around 980 MeV for all i,j = 1,2. These amplitudes are 
thus consistent with p°(770) - /o(980) mixing observed previously in the reduced transversity 



amplitudes [16(] 



The physical solutions for the relative phases LOij = vr satisfy trivially the self-consistency 
condition (2.7) that must be satisfled by co-evolution amplitudes (3.8). The four sets of physical 
solutions for transversity amplitudes Au{i), A(i{j),i, j = 1,2 can thus be identifled with the four 
sets of co-evolution amplitudes (3.8) connecting in a unique way the Kraus representation to the 
experimentally measured amplitudes. This connection validates further the view of pion creation 
process -7T~p 7r~7r"'"n as an open quantum system interacting with a quantum environment jl^. 
Information about the quantum state of the environment is accessible in part by determination 
of the probabilities pij in measurements of recoil hyperon polarization in reactions such as 
'K~p TT~K~^A^ and K~p ir^n^A^ on polarized target. 



Availability of data files from amplitude analyses. 



Output data files from computer codes for amplitude analyses are available on request to the 
author at svec@hep. physics. mcgill.ca. Code AACERNM performs amplitude analysis of mass 
dependence of the CERN data for polarized or unpolarized target. It calculates moduli and phases 
of reduced transversity amplitudes, partial wave intensities and polarizations, interference terms, 
numerous tests and auxiliary calculations. The output file for polarized target is AACERNMlp. 
The moduli and phases of reduced transversity amplitudes form input file for computer code 
HACERNMa. The code HACERNMa selects the solution n for oo and calculates first moduli 
squared and relative phases of helicity amplitudes and performs tests on their solutions before it 
calculates moduli and phases as well as real and imaginary parts of both transversity and helicity 
amplitudes in output file HACERNMan. 
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COS (CO) =-W2/W cos (CO) =+W2/W 




(7u"7u+) (GeV) 



FIG. 1: Test of consistency of the four sets of solutions Au{i), Ad{j),i,j = 1,2 with cosw = ±W2/W. The 
assumption that all four sets are physical requires that A = (wn — U21) — (wi2 — 1^22) = 0. 
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Solution set ( + +) with cos (CO) =-W2/W 
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FIG. 2: Solutions for helicity flip amplitude Li with signs of phases ++ and cosw = — W2/W. 
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Solution set ( + -) with cos (CO) =-W2/W 




0.6 0.7 0.8 0.9 1 1.1 0.6 0.7 0.8 0.9 1 1.1 



m(7r 71 ) (GeV) 

FIG. 3: Solutions for helicity flip amplitude Li with signs of phases H — and cosw = —W2/W. 



Sets (++) and ( + -) with sin(CO)=-l 
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FIG. 4: Solutions for helicity flip amplitude with sin a; = —1. The solutions are the same bor both 
sets ++ and H — of signs of phases. The helicity flip amplitude = ■^(I'S'I + and helicity non-flip 

amplitude So = -j^{\S\ + i\S\). 
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Solution set (++) with sin (CO) =-1 
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FIG. 5: Solutions for helicity flip amplitude Li with signs of phases ++ and sinw = — 1. 
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Solution set (+-) with sin (CO) =-1 
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FIG. 6: Solutions for helicity flip amplitude Li with signs of phases H — and sinw = — 1. 
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Sets (++) and ( + -) with cos(CO)=-l 
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FIG. 7: Solutions for helicity flip amplitude with cosw = — 1. The solutions are the same bor both sets 
++ and H — of signs of phases. The helicity flip amplitude = + and helicity non-flip amplitude 

So = - \S\) (Fig. 12). 
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Solution set (++) with cos (CO) =-1 
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FIG. 8: Solutions for helicity flip amplitude Li with signs of phases ++ and cosw = — 1. 
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Solution set (+-) with cos (CO) =-1 
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FIG. 9: Solutions for helicity flip amplitude Li with signs of phases H — and cosw = — 1. 
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Solution set (++) with cos (CO) =-1 
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FIG. 10: Solutions for helicity flip amplitude Ui with signs of phases ++ and cosw = — 1. 
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Solution set (+-) with cos (CO) =-1 
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FIG. 11: Solutions for helicity flip amplitude Ui with signs of phases H — and cosw = — 1. 



32 



Sets (++) and ( + -) with cos(CO)=-l 
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FIG. 12: Solutions for helicity non-flip amplitude Sq with cosw = —1. The solutions are the same bor both 
sets ++ and H — of signs of phases. The helicity non-flip amplitude = ■^(I'S'I — l^l) and helicity flip 

amplitude Si = -j={\S\ + \S\) (Fig. 7). 
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Solution set (++) with cos (CO) =-1 
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FIG. 13: Solutions for helicity non-flip amplitude Lq with signs of phases +- 1- and cos a; = — 1. 
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Solution set (+-) with cos (CO) =-1 
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FIG. 14: Solutions for helicity non-flip amplitude Lq with signs of phases H — and cos a; = — 1. 
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Solution set (++) with cos (CO) =-1 
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FIG. 15: Solutions for helicity non-flip amplitude Uq with signs of phases +- 1- and cosw = — 1. 
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Solution set (+-) with cos (CO) =-1 
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FIG. 16: Solutions for helicity non-flip amplitude Uq with signs of phases H — and cos a; = — 1. 
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Sets (++) and ( + -) with cos(CO)=-l 
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FIG. 17: Original phases of helicity flip amplitude Li for solutions with cosw = —1. Absolute phase is set 
at (i) = -TT, i = 1, 2. 
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Sets (++) and ( + -) with cos(CO)=-l 
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FIG. 18: Modified phases of helicity non-flip amplitude Sq for solutions with cosw = — 1 after assigning the 
amplitude Li the Breit-Wigner phase ^{p^) of p^{770) resonance. 
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Sets (++) and ( + -) with sin(CO)=-l 
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FIG. 19: Original phases of helicity flip amplitude Li for solutions with sinw = —1. Absolute phase is set 
at $s„(i) = -F7r/2,i = 1,2. 
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Sets (++) and ( + -) with sin(CO)=-l 
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FIG. 20: Modified phases of fielicity non-fiip amplitude Sq for solutions with sinw = —1 after assigning the 
amplitude Li the Breit-Wigner phase of p^{770) resonance. 
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